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Abstract. We have observed the high-frequency-peaked BL Lacertae object PKS 2155-304 in 2004, 2005 and 2006 with 
the CANGAROO-III imaging atmospheric Cherenkov telescope, and have detected a signal above 660 GeV at the 4.8 a level 
during the 2006 outburst period. Intranight flux variability on time scale of half an hour is observed. From this variability time 
scale, the size of the TeV gamma-ray emission region is restricted to 5 x 10'^5 cm, and the super massive black hole mass 
is estimated to be less than 1.9 x 10 SMsoiar^ where 5 is the beaming factor. The differential energy spectrum is obtained, 
and an upper limit of the extragalactic infrared background light (EBL) flux is derived under some assumption. We also fit a 
synchrotron self Compton (SSC) model to the spectral energy distribution (SED) and derive the beaming factor and magnetic 
field strength. 
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INTRODUCTION 

The first detection of TeV gamma rays from 
PKS 2155— 304 was reported by the Durham group 
during the active phase in 1997 11]. Although the 
CANGAROO-I and the CANGAROO-II observations 
in 1997, 1999, 2000, and 2001 could not find an excess 
signal from PKS 2155-304 it was confirmed 
by the H.E.S.S. group as a TeV gamma-ray source 
in 2004 [6]. In July 2006, the extreme TeV flare was 
detected by the H.E.S.S. group, and they reported very 
beautiful light curve with 1 minute timescale resolu- 
tion 17]. Some modeling for the rapid TeV flux variability 
and the spectral ener gy d istribution of PKS 2155— 304 
are reported llH |^ ITol ll 1.1 . and intensive interpretations 
of the data have been attempted. 

CANGAROO-III observations of PKS 2155-304 
have been performed from 2004 to 2006. Particularly 
in 2006, target-of-opportunity observations were made 
triggered by the H.E.S.S. report. In this paper, we 
summarize our observations and the results particularly 
focused on the data of the 2006 outburst period lfl2ll . and 
give some implications from them. 

OBSERVATIONS AND ANALYSIS 

The CANGAROO-III imaging atmospheric Cherenkov 
telescope system is operated in Woomera, South Aus- 
tralia (longitude 136°47'E, latitude 31°06'S, 160 m 
a.s.I.). Three of the four telescopes were used in these ob- 
servations. The details of the CANGAROO-III telescope 



system are described in ifTslfTilfTsll . 

Observations of PKS 2155-304 in 2004 were made 
by local trigger using three telescopes. From 2005, an 
on-line stereo trigger system became available. The de- 
tails of the CANGAROO-III trigger system are given in 
llT6llT7ll . In 2005, observations were made with only two 
telescopes due to electronics problems. Observations in 
2006 consist mostly a three-telescope configuration, but 
the exception to this is in a part of July due to a mechan- 
ical tracking problem with the third telescope. These ob- 
servations were made using wobble mode, in which the 
pointing position of each telescope was shifted in decli- 
nation by ±0.5° from the center of PKS 2155-304 al- 
ternatively every 20 minutes. Observation periods, ob- 
servation nights, the amount of observation time, the 
number of telescopes available, and the trigger mode for 
each year are summarized in TableflJ Typical trigger rate 
of the two-fold and the three-fold coincidences on July 
2006 observations were ^ 20 Hz and ^12 Hz, respec- 
tively. 

After the image cleaning and the reconstruction of the 
arrival direction, we applied the Fisher discriminant (FD) 
method to the data in order to reject numerous cosmic- 
ray background events! 18]. Following a Monte Carlo 
simulation study, the optimum FD cut values were de- 
termined, and the 9^ cuts were applied, for example at 
9^ < 0.06 deg^ for the three fold data. Details of the 
observations in 2006 and analysis method are described 
in ||12[1 . and references therein. After the data reduction 
described above, the effective live times are 17.0 hrs, 
38.6 hrs, 25.1 hrs, and 17.1 hrs in 2004, 2005, July 2006, 
and August 2006, respectively. 



TABLE 1. Summary of CANGAROO-III observations of PKS 2155-304 



Year 


Period 


Nights 


Obs. time [lirs] 


ftoftel. 


Stereo mode 


2004 


Aug.8-23 


11 


20.5 


3 


off-line 


2005 


June 6-15 


6 


46.8 


2 


on-line 




July 1-5 


5 


T 


2 


on-line 




July 29-Aug.5 


12 


T 


2 


on-line 


2006 


July 28-Aug.2 


5 


29.0 


3&2 


on-line 




Aug. 17-25 


6 


19.1 


3 


on-line 
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FIGURE 1. Long-term variations of the flux of PKS 2155- 
304 (black closed circle and bars with an arrow) compared to 
H.E.S.S. results (gray squares)[7.. ,21i1 . The flux of the H.E.S.S. 
are converted to the values above 660 GeV assuming a photon 
index of 3.3. 
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FIGURE 2. Daily light curve of PKS 2155-304 from July 
28 to August 2 in 2006 expressed by the integral flux above 
660 GeV. The dashed line indicates an average integral flux 
during this observation period. 



RESULTS 

In 2004 and 2005, no significant excess events were 
found and 2 a upper limits on the integral flux were 
obtained: 4.5 x 10-^\m-h-^ above 580 GeV for 2004 
and 6.4 x lO-'^cm-^g-i above 660 GeV for 2005. In 
July 2006, just after the largest outburst reported by 
the H.E.S.S. group|19], we detected 322 ± 67 excess 
events corresponding to 4.8a level from the direction 
of PKS2155-304||ll. The time-averaged integral flux 
above 660 GeV is calculated to be F{> 660 GeV) = 
(1.6±0.3,,«r ±0.5,,,,) X 10""cm-2s-i, which corre- 
sponds to ~ 45% of the flux observed from the Crab 
Nebula f^O*] . In the follow-up observations in August 
2006, only the flux upper limit (2 a level) was obtained: 
F{> 660 GeV) < 6.4 x lO^'^cm-^s-i (20% Crab flux), 
which means the TeV gamma-ray activity subsided in 
August. These results are summarized in Table|2] Our 
results obtained with the CANGAROO-III telescope for 
each year are compared to the H.E.S.S. results in Fig.[Tl 
where the fluxes by the H.E.S.S. are converted to the val- 
ues above 660 GeV assuming a photon index of 3.3 and 
plotted. 

Fig|2] and [3] show the nightly average integral flux 
above 660 GeV in 2006. From the light curve, it is seen 
that the averaged flux reached ~ 70% of Crab flux level 
on the night of July 30. Assuming a constant average flux 
between July 28 and August 2, a fit yields a value of 
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FIGURE 3. Daily light curve of PKS 2155-304 from August 
17 to 25 in 2006 expressed by the integral flux above 660 GeV. 
The dotted line indicates one Crab flux level above 660 GeV. 



13.9 for 4 degree of freedom, which corresponds to a 
probability of 8 x 10"^. 

We further investigate intranight variation by divid- 
ing the data into 40 minute bins for each night in 
July 2006. The light curve is shown in Fig.|4] Assum- 
ing a constant average flux for each night, test for 
flie data of July 28 and 30 gives x^/dof = 29.2/6 
and 22.1/6, respectively, which correspond to x^ prob- 
ability of 6 x 10^^ and 1 x 10^^, respectively. The in- 
tranight variations are apparent. The fractional root mean 
square variability amplitudes F,,„r are calculated to be 
0.75 ± 0.07 on July 28 and 0.58 ± 0.08 on July 30. These 
values are comparable to the variability reported by the 



TABLE 2. Summary of results for PKS2 155-304, where 2006 July corresponds to 
period from July 28 to August 2 2006 









Excess 


Significance 


Flux 


Crab flux 


Year 




[GeV] 


events 


[CT] 


[10-i2<.n,-2s-i] 


[%] 


2004 




580 


-14±28 


-0.5 


<4.5 


<10 


2005 




660 


64±53 


1.2 


<6.4 


<20 


2006 


July 


660 


322±67 


4.8 


16±3s/fl/ ±5.svs( 


45 




Aug. 


660 


-1±30 


-0.0 


<6.5 


<20 




FIGURE 4. Light curve of PKS 2155-304 expressed by the 
integral flux above 660 GeV between July 28 and August 2 
based on a 40 minute bin width. Triangles and closed circles 
indicate the flux obtained from the two-fold and three-fold data 
set, respectively. The dashed line shows an average integral flux 
during this observation period. The shaded areas indicate the 
H.E.S.S. observation periods. The time difference between the 
H.E.S.S. and CANGAROO-III sites is ~ 8 hours. 



H.E.S.S. i?!]- Doubling times are also calculated from 
Fig.|4] We obtained the shortest doubling time of 34 min- 
utes as an la upper limit. 

A photon index of the time-averaged differential 
energy spectrum of PKS 2155— 304, which was ob- 
tained between July 28 and August 2 in 2006, is 
r = —2.5 ± 0.5v,„, ± 0.7„.j,, and a flux normaliza- 
tion at 1 TeV is N{1 TeV) = (1.0 ± 0.2 ±,,„, ±,,.,.,) x 
lO-^'cm-^s-iTeV-'IT 



DISCUSSION 

Assuming that the gamma-rays from PKS 2155—304 are 
generated in the jet which is directed toward us with 
a beaming factor 8, the observed variability time scale 
fvur is related to the size R of emission region by the 
equation/? < dctvar/i^ +z) from causality, where z is the 
redshift. Using an upper limit of the shortest doubling 
time as an variability time scale, i.e. t^r — 34 minutes, 
the size of TeV gamma ray emission region is limited 
to7?<5xlO'^5 cm. The central engine of a blazar is 
believed to contain a super massive black hole. If we 
assume that the size of emission region is larger than 
the Schwarzschild radius Rsdi^ i-^- R > Rschi ^ black 



hole mass Mg is less than [c^ /2G)R. Then, substituting 
numerical values of constant, an upper limit of a super 
massive black hole mass is estimated to be Mg < 1 .9 x 
10**5Msoiar, where Msoiar is one Solar mass. On the 
other hand, assuming an internal shock model in the 
jet, where the second blob catching up with the first 
blob and create a shock wave, we estimate the distance 
D from the center of black hole to the shock region 
as D = 103(A:/10)(7/10)2[2aV(«^ - l)]Rsch- Here y 
and ay are Lorentz factors of the first and the second 
blobs(a > 1), respectively, and blobs are assumed to 
be emitted at a time interval of ^kRsch/c {k > 3). The 
size of gamma-ray emission region is D/y and observed 
variability time scale is t^ar^D / cy^ . So the black hole 
mass is expressed as 



M= l.Ox lO'' 



tvc. 



a" 



1 



2a 



Msola 



Substituting tmr = 34 minutes, and assuming a ^ I and 
k > 3, the black hole mass is estimated to be less than 

3.4 X lO^Msolar. 

As is well known, TeV photons are absorbed by 
an interaction with extragalactic infrared background 
light (EBL) through the pair production process yxev + 
ysBL^e^ + ■ Measured flux F,,, is related to the source 
flux by Fm = e^^Fs, where T is the optical depth of the 
above process. Assuming an intrinsic photon index is not 
harder than — 1.5|'23l and no absorption below 200 GeV, 
the upper limit of the optical depth t is calculated from 
the observed spectrum. Adopting the 1 a upper limit of 
the observed spectrum in July 2006, the optical depth, 
for example at 1.1/xm, is calculated to T = 2.3, and an 
upper limit on the EBL flux can be estimated to be 
45.5nWm^^sr"^ for 1.1/xm. 

On the contrary, assuming the EBL density model, 
an intrinsic spectrum can be estimated from the mea- 
sured spectrum. Choosing the EBL density model of Pri- 
mack et al. ll23ll . an intrinsic spectrum of dN/dE = (2.7 ± 



-(1.8±0.6) 



.-2,-1- 



TeV ^ is de- 



0.6) X 10-"(£/ireV) 
rived from the observed spectrum. When applying the 
model by Stecker et al. 12411 to the observed spectrum, the 
deabsorbed spectrum seems to be too hard. 

Using the intrinsic energy spectrum derived above, 
the simple SSC model by Kino et al. 1250 is applied 
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FIGURE 5. Spectral energy distribution of PKS 2 1 55-304 in 
2006. Red butterfly area indicates the deabsorbed average flux 
above 660 GeV observed with the CANGAROO-III telescope 
between July 28 and August 2, and light blue shows the results 
by H.E.S.S. on July 27|7]. Green plots indicate the data of XTE 
and UVOT on Swift] 8] and Narrabri on July 30. Deep green and 
black indicate Swift data on Aug.2 and RXTE data on July 30, 
respectively. 

to the spectral energy distribution. Figure |5] shows the 
spectral energy distributions of multiwavelength obser- 
vations. Adopting 5 x 10'^5 cm as the radius of the emis- 
sion region based on our results, fitting was done to our 
average spectral energy distribution between July 28 and 
August 2 2006 and X-ray data taken by Swift on July 
30, which are not perfectly simultaneous observations. 
We found that the beaming factor 5 ^ 60 and magnetic 
field B ~ 2.5 G explain the observed spectral energy dis- 
tribution well. This beaming factor is much greater than 
the typical value for blazars and is consistent with values 
obtained by other (§, . However the magnetic field is 
stronger compared to those they derived. We also try to 
fit the H.E.S.S. data and we found the smaller size of 
emission region and higher electron density are required, 
and beaming factor 5 ^ 80 and magnetic field B 5.0G 
are obtained. However, it is still poorly-fitting, and we 
need to consider more complicated and realistic model. 

CONCLUSION 

We observed the nearby HBL PKS 2155-304 in 2004, 
2005, and 2006 with the CANGAROO-III imaging at- 
mospheric Cerenkov telescope. During the active phase 
in July 2006, we detected a signal at 4.8 (7 level 
45 % Crab flux) above 660 GeV. The intranight time vari- 
ations of the flux were observed, and obtained shortest 
doubling time was 34 minutes. The size of TeV gamma- 
ray emission region is limited to 5 x 10'^5 cm from the 
variability time scale, and a super massive black hole 
mass is restricted to be less than 1 .9 x 10*^5Msoiar, where 



5 is a beaming factor. From the differential energy spec- 
trum, an upper limit of the extragalactic infrared back- 
ground light (EBL) flux is derived under the model as- 
sumption. The synchrotron self Compton (SSC) model 
fitting suggests a larger beaming factor compared to typ- 
ical blazars and a rather larger magnetic strength. 
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